Biodiversity Medicine: New Horizon and New Opportunity for Cancer.

SUMMARY: Accessibility to standard of care remains a challenge to patients in low- and middle-income countries (LMIC), hampering efforts to alleviate the burden of cancer and to improve overall health outcomes. In response to this pressing global health care issue, we propose here a new strategy to create affordable, easily accessible, and effective therapeutic solutions to address this inequity in cancer treatment: the use of science-based biodiversity medicine as an alternative to modern drug therapy, in which we will leverage and combine high-throughput omics technologies with artificial intelligence, to study local biodiversity, their potential anticancer properties, and short- and long-term clinical response and outcomes.

Establishment and characterization of a patient-derived solitary fibrous tumor/hemangiopericytoma cell line model.

Solitary fibrous tumor/Hemangiopericytoma (SFT/HPC) is a rare subtype of soft tissue sarcoma harboring NAB2-STAT6 gene fusions. Mechanistic studies and therapeutic development on SFT/HPC are impeded by scarcity and lack of system models. In this study, we established and characterized a novel SFT/HPC patient-derived cell line (PDC), SFT-S1, and screened for potential drug candidates that could be repurposed for the treatment of SFT/HPC. Immunohistochemistry profiles of the PDC was consistent with the patient's tumor sample (CD99+/CD34+/desmin-). RNA sequencing, followed by Sanger sequencing confirmed the pathognomonic NAB2exon3-STAT6exon18 fusion in both the PDC and the original tumor. Transcriptomic data showed strong enrichment for oncogenic pathways (epithelial-mesenchymal transition, FGF, EGR1 and TGFß signaling pathways) in the tumor. Whole genome sequencing identified potentially pathogenic somatic variants such as MAGEA10 and ABCA2. Among a panel of 14 targeted agents screened, dasatinib was identified to be the most potent small molecule inhibitor against the PDC (IC50, 473 nM), followed by osimertinib (IC50, 730 nM) and sunitinib (IC50, 1765 nM). Methylation profiling of the tumor suggests that this specific variant of SFT/HPC could lead to genome-wide hypomethylation. In conclusion, we established a novel PDC model of SFT/HPC with comprehensive characterization of its genomic, epigenomic and transcriptomic landscape, which can facilitate future preclinical studies of SFT/HPC, such as in vitro drug screening and in vivo drug testing.

Integrative multiomics enhancer activity profiling identifies therapeutic vulnerabilities in cholangiocarcinoma of different etiologies.

OBJECTIVES: Cholangiocarcinoma (CCA) is a heterogeneous malignancy with high mortality and dismal prognosis, and an urgent clinical need for new therapies. Knowledge of the CCA epigenome is largely limited to aberrant DNA methylation. Dysregulation of enhancer activities has been identified to affect carcinogenesis and leveraged for new therapies but is uninvestigated in CCA. Our aim is to identify potential therapeutic targets in different subtypes of CCA through enhancer profiling. DESIGN: Integrative multiomics enhancer activity profiling of diverse CCA was performed. A panel of diverse CCA cell lines, patient-derived and cell line-derived xenografts were used to study identified enriched pathways and vulnerabilities. NanoString, multiplex immunohistochemistry staining and single-cell spatial transcriptomics were used to explore the immunogenicity of diverse CCA. RESULTS: We identified three distinct groups, associated with different etiologies and unique pathways. Drug inhibitors of identified pathways reduced tumour growth in in vitro and in vivo models. The first group (ESTRO), with mostly fluke-positive CCAs, displayed activation in estrogen signalling and were sensitive to MTOR inhibitors. Another group (OXPHO), with mostly BAP1 and IDH-mutant CCAs, displayed activated oxidative phosphorylation pathways, and were sensitive to oxidative phosphorylation inhibitors. Immune-related pathways were activated in the final group (IMMUN), made up of an immunogenic CCA subtype and CCA with aristolochic acid (AA) mutational signatures. Intratumour differences in AA mutation load were correlated to intratumour variation of different immune cell populations. CONCLUSION: Our study elucidates the mechanisms underlying enhancer dysregulation and deepens understanding of different tumourigenesis processes in distinct CCA subtypes, with potential significant therapeutics and clinical benefits.

Classification of HER2-negative breast cancers by ERBB2 copy number alteration status reveals molecular differences associated with chromosome 17 gene aberrations.

Background: Recently, HER2-negative breast cancers have been reclassified by protein expression into 'HER2-low' and 'HER2-zero' subgroups, but the consideration of HER2-low breast cancer as a distinct biological subtype with differing prognoses remains controversial. By contrast, non-neutral ERBB2 copy number alteration (CNA) status is associated with inferior survival outcomes compared to ERBB2 CNA-neutral breast cancer, providing an alternative approach to classification. Methods: Here, we investigated the molecular landscape of non-metastatic HER2-negative BCs in relation to ERBB2 CNA status to elucidate biological differences. Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) TCGA-BRCA datasets (n = 1875) were analyzed. Results: Nearly two-fifths of the cohort harbored ERBB2 CNAs (39.4%), which were significantly enriched within hormone receptor-negative (56.1%) than within hormone receptor-positive BCs (35.5%; p < 0.0001). Globally, CNAs across the genome were significantly higher in ERBB2 non-neutral compared to neutral cohorts (p < 0.0001). Notably, genetic aberrations on chromosome 17 - BRCA1, NF1, TP53, MAP2K4, and NCOR1 - were widespread in the ERBB2 non-neutral cases. While chromosome 17q arm-level alterations were largely in tandem with ERBB2 CNA status, arm-level loss in chromosome 17p was prevalent regardless of ERBB2 gain, amplification, or loss. Differential gene expression analysis demonstrated that pathways involved in the cell cycle, proteasome, and DNA replication were upregulated in ERBB2 non-neutral cases. Conclusion: Classification of HER2-negative BCs according to ERBB2 CNA status reveals differences in the genomic landscape. The implications of concurrent aberrations in other genes on chromosome 17 merit further research in ERBB2 non-neutral BCs.

Multi-omic profiling and real time ex vivo modelling of imatinib-resistant dermatofibrosarcoma protuberans with fibrosarcomatous transformation.

Dermatofibrosarcoma protuberans (DFSP) is a rare and indolent cutaneous sarcoma, with the risk of aggressive fibro-sarcomatous transformation. Limited effective options are available for un-resectable or metastatic DFSP beyond targeting the oncogenic PDGF pathway with imatinib therapy. We established a patient-derived xenograft (PDX) and cell line model (designated MDFSP-S1) of imatinib-resistant DFSP with fibro-sarcomatous transformation. Whole genome sequencing identified high-level amplification at chromosomes 17 and 22, whilst homozygous deep deletion was demonstrated at chromosome 9 (CDKN2A, CDKN2B, MTAP). RNA sequencing followed by Sanger sequencing confirmed the pathognomonic COL1A1-PDGFB t (17;22) rearrangement in the original tumour, PDX and cell line model. Immunohistochemistry profiles of the PDX model were consistent with the patient's tumour sample (CD34 + /MIB1 + /SOX10- ). Gene set enrichment analysis highlighted top-scoring Hallmark gene sets in several oncogenic signalling pathways, including potentially targetable MTORC1 signalling and angiogenesis pathways. Antiangiogenic agents (sunitinib, regorafenib, pazopanib, axitinib) and the third-generation irreversible epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor osimertinib exhibited modest anti-proliferative activity in the cell line, with IC50 values between 1 and 10 µM at 72 h. No significant activity was observed with imatinib, palbociclib, everolimus, olaparib, gefitinib and erlotinib (IC50 all > 10 µM). In conclusion, we established MDFSP-S1, a new PDX and cell line model of imatinib-resistant DFSP with fibro-sarcomatous transformation.

Giant juvenile fibroadenomas with and without prominent pseudoangiomatous stromal hyperplasia (PASH)-like change: clinicopathological and molecular characteristics.

AIMS: Juvenile fibroadenomas (JFA) are biphasic fibroepithelial lesions (FEL) usually occurring in adolescent female patients. Giant (G) JFA, like other FEL, may exhibit prominent pseudoangiomatous stromal hyperplasia (PASH)-like change. We sought to determine clinicopathological and molecular characteristics of GJFA with and without PASH. METHODS AND RESULTS: Archives were searched for cases of GJFA (1985-2020). All were stained for androgen receptor (AR), beta-catenin, CD34 and progesterone receptor (PR). Cases were sequenced using a custom 16-gene panel - MED12 (exons 1 and 2), TERT promoter (-124C>T and -146Ctable>T), SETD2, KMT2D, RARA (exons 5-9), FLNA, NF1, PIK3CA (exons 10, 11 and 21), EGFR, RB1, BCOR, TP53, PTEN, ERBB4, IGF1R and MAP3K1. Twenty-seven GJFA from 21 female patients aged 10.1-25.2 years were identified. Size ranged from 5.2 to 21 cm. Two patients had multiple, bilateral and later recurrent GJFA. Thirteen (48%) cases showed prominent PASH-like stroma. All were positive for stromal CD34, negative for AR and beta-catenin and one case showed focal PR expression. Sequencing showed MAP3K1 and SETD2 mutations in 17 samples, with KMT2D, TP53 and BCOR aberrations in 10 (45%), 10 (45%) and seven (32%) cases, respectively. Tumours with a PASH-like pattern had higher prevalence of SETD2 (P = 0.004) and TP53 (P = 0.029) mutations, while those without PASH had more RB1 mutations (P = 0.043). MED12 mutation was identified in one case. TERT promoter mutation was observed in four (18%), including two recurrences. CONCLUSIONS: Gene mutations along more advanced phases of the proposed FEL pathogenetic pathway in GJFA are unusual, and suggest a mechanism for more aggressive growth in these tumours.

Spatial transcriptomics reveal topological immune landscapes of Asian head and neck angiosarcoma.

Angiosarcomas are rare malignant tumors of the endothelium, arising commonly from the head and neck region (AS-HN) and recently associated with ultraviolet (UV) exposure and human herpesvirus-7 infection. We examined 81 cases of angiosarcomas, including 47 cases of AS-HN, integrating information from whole genome sequencing, gene expression profiling and spatial transcriptomics (10X Visium). In the AS-HN cohort, we observed recurrent somatic mutations in CSMD3 (18%), LRP1B (18%), MUC16 (18%), POT1 (16%) and TP53 (16%). UV-positive AS-HN harbored significantly higher tumor mutation burden than UV-negative cases (p = 0.0294). NanoString profiling identified three clusters with distinct tumor inflammation signature scores (p < 0.001). Spatial transcriptomics revealed topological profiles of the tumor microenvironment, identifying dominant but tumor-excluded inflammatory signals in immune-hot cases and immune foci even in otherwise immune-cold cases. In conclusion, spatial transcriptomics reveal the tumor immune landscape of angiosarcoma, and in combination with multi-omic information, may improve implementation of treatment strategies.

A comprehensive next generation sequencing tissue assay for Asian-prevalent cancers-Analytical validation and performance evaluation with clinical samples.

Introduction: A well-validated diagnostic assay with curated biomarkers complements clinicopathological factors to facilitate early diagnosis and ensure timely treatment delivery. This study focuses on an Asian-centric cancer diagnostic assay designed and thoroughly validated against commercially available standard references and a cohort of over 200 clinical specimens spanning 12 diverse Asian-centric cancer types. Methods: The assay uses hybrid-capture probes capable of profiling DNA aberrations from 572 cancer-related genes and 91 RNA fusion partners. The panel can detect clinically-tractable biomarkers such as microsatellite instability (MSI) and tumor mutation burden (TMB). Results: Analytical evaluation demonstrated 100% specificity and 99.9% sensitivity within a ≥5% VAF limit of detection (LoD) for SNV/Indels. RNA-based fusion features an LoD of ≥5 copies per nanogram input when evaluated against commercial references. Excellent linearity and concordance were observed when benchmarking against orthogonal methods in identifying MSI status, TMB scores and RNA fusions. Actionable genetic alterations were identified in 65% of the clinical samples. Conclusion: These results demonstrate a molecular diagnostic assay that accurately detects genomic alterations and complex biomarkers. The data also supports an excellent performance of this assay for making critical diagnoses and well-informed therapeutic decisions in Asian prevalent cancers.

The Multi-Dimensional Biomarker Landscape in Cancer Immunotherapy.

The field of immuno-oncology is now at the forefront of cancer care and is rapidly evolving. The immune checkpoint blockade has been demonstrated to restore antitumor responses in several cancer types. However, durable responses can be observed only in a subset of patients, highlighting the importance of investigating the tumor microenvironment (TME) and cellular heterogeneity to define the phenotypes that contribute to resistance as opposed to those that confer susceptibility to immune surveillance and immunotherapy. In this review, we summarize how some of the most widely used conventional technologies and biomarkers may be useful for the purpose of predicting immunotherapy outcomes in patients, and discuss their shortcomings. We also provide an overview of how emerging single-cell spatial omics may be applied to further advance our understanding of the interactions within the TME, and how these technologies help to deliver important new insights into biomarker discovery to improve the prediction of patient response.

Therapeutic and immunomodulatory potential of pazopanib in malignant phyllodes tumor.

Malignant phyllodes tumors (PT) are rare aggressive fibroepithelial neoplasms with high metastatic potential and lack effective therapy. We established a patient-derived xenograft (PDX) and cell line model (designated MPT-S1) of malignant PT which demonstrated clinical response to pazopanib. Whole exome sequencing identified somatic mutations in TP53, RB1, MED12, and KMT2D. Immunohistochemistry and genomic profiles of the tumor, PDX and cell line were concordant. In keeping with clinical observation, pazopanib reduced cell viability in a dose-dependent manner and evoked apoptosis, and led to significant abrogation of in vivo tumor growth. Whole transcriptomic analysis revealed that pazopanib decreased expression of genes involved in oncogenic and apoptosis signaling. We also observed decreased expression of ENPP1, with known roles in cancer invasion and metastasis, as well as STING pathway upregulation. Accordingly, pazopanib induced micronuclei formation, and evoked phospho-TBK1 and PD-L1 expression. In an additional cohort of malignant PT (n = 14), six (42.9%) showed comparable or higher levels of ENPP1 relative to MPT-S1, highlighting its potential role as a therapeutic target. In conclusion, we established MPT-S1, a new PDX and cell line model, and provided evidence for the clinical efficacy of pazopanib in malignant PT.

HER2 expression, copy number variation and survival outcomes in HER2-low non-metastatic breast cancer: an international multicentre cohort study and TCGA-METABRIC analysis.

BACKGROUND: HER2-low breast cancer (BC) is currently an area of active interest. This study evaluated the impact of low expression of HER2 on survival outcomes in HER2-negative non-metastatic breast cancer (BC). METHODS: Patients with HER2-negative non-metastatic BC from 6 centres within the Asian Breast Cancer Cooperative Group (ABCCG) (n = 28,280) were analysed. HER2-low was defined as immunohistochemistry (IHC) 1+ or 2+ and in situ hybridization non-amplified (ISH-) and HER2-zero as IHC 0. Relapse-free survival (RFS) and overall survival (OS) by hormone receptor status and HER2 IHC 0, 1+ and 2+ ISH- status were the main outcomes. A combined TCGA-BRCA and METABRIC cohort (n = 1967) was also analysed to explore the association between HER2 expression, ERBB2 copy number variation (CNV) status and RFS. RESULTS: ABCCG cohort median follow-up was 6.6 years; there were 12,260 (43.4%) HER2-low BC and 16,020 (56.6%) HER2-zero BC. The outcomes were better in HER2-low BC than in HER2-zero BC (RFS: centre-adjusted hazard ratio (HR) 0.88, 95% CI 0.82-0.93, P < 0.001; OS: centre-adjusted HR 0.82, 95% CI 0.76-0.89, P < 0.001). On multivariable analysis, HER2-low status was prognostic (RFS: HR 0.90, 95% CI 0.85-0.96, P = 0.002; OS: HR 0.86, 95% CI 0.79-0.93, P < 0.001). These differences remained significant in hormone receptor-positive tumours and for OS in hormone receptor-negative tumours. Superior outcomes were observed for HER2 IHC1+ BC versus HER2-zero BC (RFS: HR 0.89, 95% CI 0.83-0.96, P = 0.001; OS: HR 0.85, 95% CI 0.78-0.93, P = 0.001). No significant differences were seen between HER2 IHC2+ ISH- and HER2-zero BCs. In the TCGA-BRCA and METABRIC cohorts, ERBB2 CNV status was an independent RFS prognostic factor (neutral versus non-neutral HR 0.71, 95% CI 0.59-0.86, P < 0.001); no differences in RFS by ERBB2 mRNA expression levels were found. CONCLUSIONS: HER2-low BC had a superior prognosis compared to HER2-zero BC in the non-metastatic setting, though absolute differences were modest and driven by HER2 IHC 1+ BC. ERBB2 CNV merits further investigation in HER2-negative BC.

Ligand-mediated PAI-1 inhibition in a mouse model of peritoneal carcinomatosis.

Peritoneal carcinomatosis (PC) present a ubiquitous clinical conundrum in all intra-abdominal malignancies. Via functional and transcriptomic experiments of ascites-treated PC cells, we identify STAT3 as a key signaling pathway. Integrative analysis of publicly available databases and correlation with clinical cohorts (n = 7,359) reveal putative clinically significant activating ligands of STAT3 signaling. We further validate a 3-biomarker prognostic panel in ascites independent of clinical covariates in a prospective study (n = 149). Via single-cell sequencing experiments, we uncover that PAI-1, a key component of the prognostic biomarker panel, is largely secreted by fibroblasts and mesothelial cells. Molecular stratification of ascites using PAI-1 levels and STAT3 activation in ascites-treated cells highlight a therapeutic opportunity based on a phenomenon of paracrine addiction. These results are recapitulated in patient-derived ascites-dependent xenografts. Here, we demonstrate therapeutic proof of concept of direct ligand inhibition of a prognostic target within an enclosed biological space.

Rare Occurrence of Aristolochic Acid Mutational Signatures in Oro-Gastrointestinal Tract Cancers.

BACKGROUND: Aristolochic acids (AAs) are potent mutagens commonly found in herbal plant-based remedies widely used throughout Asian countries. PATIENTS AND METHODS: To understand whether AA is involved in the tumorigenesis of the oro-gastrointestinal tract, we used whole-exome sequencing to profile 54 cases of four distinct types of oro-gastrointestinal tract cancer (OGITC) from Taiwan. RESULTS: A diverse landscape of mutational signatures including those from DNA mismatch repair and reactive oxygen species was observed. APOBEC mutational signatures were observed in 60% of oral squamous cell carcinomas. Only one sample harbored AA mutational signatures, contradictory to prior reports of cancers from Taiwan. The metabolism of AA in the liver and urinary tract, transient exposure time, and high cell turnover rates at OGITC sites may explain our findings. CONCLUSION: AA signatures in OGITCs are rare and unlikely to be a major contributing factor in oro-gastrointestinal tract tumorigenesis.

Multiregion sequencing of sarcomatoid renal cell carcinoma arising from autosomal dominant polycystic kidney disease.

BACKGROUND: Autosomal dominant polycystic kidney disease (ADPKD) is an inherited cystic kidney disease associated with a spectrum of various renal and extrarenal manifestations, including increased risk of kidney cancers. Here, we present the initial molecular description of sarcomatoid renal cell carcinoma (sRCC) arising in the setting of ADPKD. METHODS: Multiregion whole-exome sequencing and whole transcriptomic sequencing were used to examine intratumoral molecular heterogeneity among histologically-distinct spindle (sarcomatoid), epithelioid, or biphasic compartments within the tumor and compared with the non-malignant ADPKD component. RESULTS: Spindle and biphasic components harbored several overlapping driver gene mutations, but do not share any with the epithelioid component. Mutations in ATM, CTNNB1, and NF2 were present only in the biphasic and spindle components, while mutations in BID, FLT3, ARID1B, and SMARCA2 were present only in the epithelioid component. We observed dichotomous evolutionary pathways in the development of epithelioid and spindle compartments, involving early mutations in TP53 and ATM/CTNNB1/NF2 respectively. Wnt, PI3K-mTOR, and MAPK signaling pathways, known key mechanisms involved in ADPKD development, featured prominently in the sarcomatoid component. CONCLUSION: This highlights that common pro-oncogenic signals are present between ADPKD and sRCC providing insights into their shared pathobiology.

Genetic differences between benign phyllodes tumors and fibroadenomas revealed through targeted next generation sequencing.

Breast fibroepithelial lesions are biphasic tumors which comprise the common benign fibroadenomas (FAs) and the rarer phyllodes tumors (PTs). This study analyzed 262 (42%) conventional FAs, 45 (7%) cellular FAs, and 321 (51%) benign PTs contributed by the International Fibroepithelial Consortium, using a previously curated 16 gene panel. Benign PTs were found to possess a higher number of mutations, and higher rates of cancer driver gene alterations than both groups of FAs, in particular MED12, TERT promoter, RARA, FLNA, SETD2, RB1, and EGFR. Cases with MED12 mutations were also more likely to have TERT promoter, RARA, SETD2, and EGFR. There were no significant differences detected between conventional FAs and cellular FAs, except for PIK3CA and MAP3K1. TERT promoter alterations were most optimal in discriminating between FAs and benign PTs. Our study affirms the role of sequencing and key mutations that may assist in refining diagnoses of these lesions.

Genomic Landscapes of Acral Melanomas in East Asia.

BACKGROUND/AIM: Acral melanomas (AM) represent a rare subgroup of melanomas with poor clinical outcomes and are enriched in Asian populations. Recent advances in next generation sequencing have provided opportunities to apply precision medicine to AM. PATIENTS AND METHODS: Here, we present a series of 13 patients with melanomas from Taiwan and Singapore, including 8 patients with AM profiled using whole exome sequencing and summarize the recent studies on the genomic landscape of AM. RESULTS: We identified mutually exclusive mutations in BRAF, NRAS, HRAS, NF1 and KIT in 6 AM cases. In addition, recurrent copy number gains in CCND1 and CDK4, as well as recurrent deletions in CDKN2A/CDKN2B, ATM and RAD51 were observed, supporting the potential use of CDK4/6 or PARP inhibitors in the treatment of these patients. CONCLUSION: The genomic landscape of AM provides an important resource for applying novel targeted therapies in this rare disease.

Bizarre giant cells in human angiosarcoma exhibit chemoresistance and contribute to poor survival outcomes.

Giant cells (GC) are a poorly understood subset of tumor cells that have been increasingly recognized as a potential contributor to tumor heterogeneity and treatment resistance. We aimed to characterize the biological and clinical significance of GC in angiosarcoma, an aggressive rare cancer of endothelial origin. Archival angiosarcoma samples were examined for the presence of GC and compared with clinicopathological as well as NanoString gene expression data. GC were examined in angiosarcoma cell lines MOLAS and ISOHAS using conventional and electron microscopy, single cell whole genome profiling, and other assays. In the cell lines, GC represented a rare population of mitotically active, non-senescent CD31+ cells, and shared similar genomic profiles with regular-sized cells, consistent with a malignant endothelial phenotype. GC remained viable and persisted in culture following exposure to paclitaxel and doxorubicin. In patient samples, GC were present in 24 of 58 (41.4%) cases. GC was correlated with poorer responses to chemotherapy (25.0% vs 73.3%, P = 0.0213) and independently contributed to worse overall survival outcomes (hazard ratio 2.20, 95% confidence interval 1.17-4.15, P = 0.0142). NanoString profiling revealed overexpression of genes, including COL11A1, STC1, and ERO1A, accompanied by upregulation of immune-related metabolic stress and metastasis/matrix remodeling pathways in GC-containing tumors. In conclusion, GC may contribute to chemoresistance and poor prognosis in angiosarcoma.

PSiTE: a Phylogeny guided Simulator for Tumor Evolution.

SUMMARY: Simulating realistic clonal dynamics of tumors is an important topic in cancer genomics. Here, we present Phylogeny guided Simulator for Tumor Evolution, a tool that can simulate different types of tumor samples including single sector, multi-sector bulk tumor as well as single-cell tumor data under a wide range of evolutionary trajectories. Phylogeny guided Simulator for Tumor Evolution provides an efficient tool for understanding clonal evolution of cancer. AVAILABILITY AND IMPLEMENTATION: PSiTE is implemented in Python and is available at https://github.com/hchyang/PSiTE. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.